Search for conformal invariance in compressible two-dimensional turbulence 
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We present a search for conformal invariance in vorticity isolines of two-dimensional compressible 
turbulence. The vorticity is measured by tracking the motion of particles that float at the surface 
of a turbulent tank of water. The three-dimensional turbulence in the tank has a Taylor microscale 
Re\ — 160. The conformal invariance theory being tested here is related to the behavior of equilib- 
rium systems near a critical point. This theory is associated with the work of Lowner, Schramm and 
others and is usually referred to as Schramm-Lowner Evolution (SLE). The system was exposed to 
several tests of SLE. The results of these tests suggest that zero- vorticity isolines exhibit noticeable 
departures from this type of conformal invariance. 
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I. INTRODUCTION 

Strong turbulence is fiercely difficult to understand be- 
cause it is dominated by nonlinear effects, and because 
many degrees of freedom of fluid flow are excited 
Some of its main features can be understood from di- 
mensional arguments based on an eddy cascade theory. 
There is a range of eddy sizes over which the system 
exhibits approximate self-similarity or scale invariance. 
Most efforts have been focused on this self-similar range, 
often called the inertial range. It can span length scales 
ranging from many meters down to eddy sizes of tens of 
microns in the atmosphere, in the ocean, and in large 
wind tunnels. 

Though most endeavors are focused on understand- 
ing three-dimensional turbulence, two dimensional (2D) 
flows are important from a practical and fundamental 
point of view. The depth of oceans (L) and the thickness 
of the atmosphere is very small compared to the earth's 
radius R. Large-scale velocity variations R >> L are 
properly viewed as two-dimensional. 

Two dimensional turbulence displays striking differ- 
ences from its 3D counterpart. In 2D, smaller eddies 
combine to form bigger ones while the reverse happens 
in three dimensions. This so-called "inverse cascade" in 
2D turbulence characterizes hurricane growth in the tro- 
posphere Q. 

The theory of two-dimensional (2D) turbulence brings 
new complications and simplifications at the same time. 
Between the large eddies in 2D turbulence are thin re- 
gions where the vorticity of the flow is very large, even 
though these regions contain only a small fraction of the 
turbulent energy. In two dimensions, the vorticity is, of 
course, perpendicular to the plane of the 2D flow and is 
hence a scalar. In 3D, the vorticity is amplified by ve- 
locity gradients, whereas in 2D, its mean square vorticity 
is a constant, viscous damping aside. As for this damp- 
ing, it occurs at small scales and is almost absent in the 



self-similar range of interest here. 

It is obviously important if 2D turbulence should turn 
out to exhibit invariance features that go beyond self- 
similarity. Recent theoretical and numerical work sug- 
gests that this is so. This evidence comes from the stud 
of contours of zero vorticity in an incompressible flow 
(to be called BBCF) . In that important paper, the au- 
thors focus on the geometry of (contorted) paths through 
the fluid where the vorticity ui is zero at each instant 
of time r. Their simulations reveal a new type of con- 
formal invariance. The conformal invariance discussed 
here , and by BBCF, is unrelated to the usual conformal 
mapping technique applied to electric potential functions. 
Rather, it is about the growth of a random curve where 
each incremental length is produced by a conformal map 
characterized by a Brownian function. 

A path in the x, y plane, written as z — (x, iy), is con- 
formally invariant (in this sense Q), if there is a function 
g{z) = (u, iv) that can map the path back to the real axis 
in the u, v plane while preserving all the angles (See Fig. 
®- 

The present work is an experimental study of the 
contours of constant vorticity in a compressible flow. 
We know of no published experiments for incompress- 
ible systems. The present laboratory observations dis- 
play approximate conformal invariance for the contours 
of zero vorticity measured at hundreds of instants of time 
r. Contours of nonzero vorticity paths were also exam- 
ined. Larger deviations from conformal invariance are 
observed. 

It was noted by BBCF that with present technology, it 
is technically not possible to search for conformal invari- 
ance in incompressible 2D flow experiments. However, in 
the compressible flow experiments to be described here, 
it was possible to accumulate data over a sufficiently wide 
parameter range to make such a test. In this work, we use 
an overhead fast camera to track the motion of particles 
that float on a turbulent tank of water. These particles 
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have a density that is a fourth of the density of water, 
so their motion is confined to the surface of the underly- 
ing turbulent flow, which is, of course three-dimensional. 
The flow of the turbulent water underneath the floaters 
is incompressible, assuring that the two-dimensional di- 
vergence of the velocity of the floaters is not zero. The 
particles are small enough to be almost inertia-free. Thus 
they sample the velocity of the flow v(x, y, z) in the sur- 
face plane z = 0, as discussed below. The present ex- 
periments are performed at a moderately high Reynolds 
number where the inertial range is appreciably large. 

This system of floaters is very different from conven- 
tional two-dimensional turbulence and also from those 
studied by BBCF. It is not merely that the surface on 
which the floaters move is rippled (their amplitude is 
small 0), it is that the floaters do not form a separate 
system; they can exchange energy with those water par- 
ticles beneath them. In principle, at least, they can take 
and return their kinetic energy to the underlying fluid 
on all spatial scales. Thus there is no reason to expect 
an energy cascade which implies dissipation only at small 
scales or at the boundaries of the container. Nevertheless, 
the squared velocity difference (5v(r) 2 ) = D2(r) between 
pairs of points separated by distances r, closely conform 
to that of 3D turbulence, D2(r) oc r 2 / 3 . This scaling is 
seen in experiments as well as simulations [5|. 

The present study hinges on the measurement of a ran- 
dom variable U(t), yet to be defined, whose average mean 
square must be Brownian in character. In that case 



W(t)-u(o)f) = K t a . 



(1) 



The exponent a and the value of k are measured in the 
experiments discussed here. 

The parameter t is a dimcnsionless length and not 
time, and the exponent a must be unity, as in Brown- 
ian motion, a requirement that must be met if confor- 
mal invariance is realized. The dimensionless "diffusiv- 
ity" k is a very important parameter in the theory be- 
ing tested here. One may think of k as the dimension- 
less diffusivity, but only if a = 1. For a self- avoiding 
random walk k =8/3, and for critical percolation k = 
6 Q. This last value is deduced in the simulations of 
BBCF. In the present experiments, the value of k for the 
zero-vorticity contours was extracted from measurements 
made at many instants of time. At each instant of time r, 
there are many constant vorticity lines. The parameter 
t increases along each line. 



II. THE SEARCH FOR U(t) 




FIG. 1: A typical SLE trace with n = 2 is shown in Fig. 
1. Writing gt{z) = u + iv,z — x + iy, the horizontal and 
vertical axes in this figure are x and iy respectively. The 
solid irregular line, which corresponds to a particular value 
the parameter t, separates a pair of regions where g(z) is an- 
alytic. It is called a trace. The traces are self-similar and also 
self-avoiding. In the experiments to be discussed below, the 
traces were measured at many of instants of (dimensionless) 
time r. (With permission from T. Kennedy. See T. Kennedy, 
|http://www.math.arizona.edu/^t gk/rtg_20 1f /sle2.0.pdf| 
2011 for original graph.) 



2D curve, which is being mapped by gt{z). This depen- 
dence on t is what is usually referred to as the Lowner 
differential equation, 



dgt(z) 
dt 



g t (z)-U(t) 



(2) 



It was O. Schramm who first discovered that for a con- 
formally invariant random curve in a 2D plane, U(t) is a 
one-dimensional Brownian motion obeying Eq. [T] [6| . In 
this case, the random curve is referred to as a Schramm- 
Lowner Evolution (SLE) trace. For a more formal and 
complete discussion of SLE, see Q. 

One of the most effective ways to identify this type of 
conformal invariance is to measure U (t) , defined by the 
above differential equation, and see if its mean square 
average obeys Eq. [T] In the next two sections, the ex- 
periment and the procedure to calculate ((SU(t)) 2 ) = 
((U(t) - U(0)) 2 ) is described. 

A typical SLE trace is shown in Fig. [TJ Such traces are 
both self-similar and self-avoiding. After the conformal 
transformation, the grid in the u, v plane is rectangular; 
the transformation under gt{z) is conformal, as all the 
angles are preserved. 



The analysis of the experimental data, to be discussed 
below, requires that the above-defined g{z) be uniquely 
determined by the function U(t) (usually referred to as 
the driving function) , which is related to the experimen- 
tal observations. Thus, g(z) also depends on t and hence 
will be written as gt{z). Roughly speaking, the dimen- 
sionless parameter t is proportional to the "length" of the 



III. EXPERIMENTAL 

The fm x I m tank is filled with water to a height of 
30 cm. The tank is large compared to the camera's field 
of view. The turbulence is generated by a large pump 
connected to a network of rotating jets in a plane 10 
cm above the tank floor. See Fig. [5] for a schematic of 
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FIG. 2: Schematic of the top-view (top panel) and side-view 
(bottom panel) of the experimental setup. 36 rotating capped 
jets are placed horizontally on the tank floor (shown as ran- 
domly oriented Z-shaped patterns) that pump water into the 
tank recirculated by an 8 hp pump. The region in the lateral 
center of the tank and at the surface (z=0) is illuminated by 
a laser-sheet. A high-speed digital camera suspended verti- 
cally above this central region captures images of the light 
scattered by buoyant particles (50 (im hollowglass spheres of 
specific gravity 0.25). 



the experimental setup. The arrangement creates uni- 
form turbulence in the center of the tank and also moves 
the source of turbulent injection far from the fluid sur- 
face where the measurements are made |5|. With this 
scheme, surface waves, which cannot be avoided, do not 
exceed an amplitude of ~ 1 mm [5j|. It is necessary that 
the surface of the tank be freshly cleaned before each set 
of measurements. Otherwise, amphiphiles form a contin- 
uous layer on the surface and prevents the floaters from 
moving freely under the action of the turbulence |5|. 

The hydrophilic particles chosen here are subject to 
capillary forces which are very small compared to forces 
coming from the turbulence, and do not affect the results 
as they do in 0, [1| . The non-inertial character of the par- 
ticles is minimal because the Stokes number St is small: 



TABLE I: Turbulent parameters measured at the surface. 
Measurements are made at several values of Re\ with an av- 
erage Rex ~ 160. The parameters listed are averages, with 
deviations less than 10%. 



Parameter 



Taylor microscale A (cm) 

Taylor Re\ 

Integral scale lo (cm) 

Large Eddy Turnover 
Time (LETT) to (s) 

Dissipation rate 
Ediss (cm 2 /s 3 ) 

Kolmogorov scale 77 (cm) 

RMS velocity iv ms (cm/s) 

Compressibility C 



Symbol used 
in text 



Rex = 



lo = J dr ", r „„r^V 



lo 



TO = 



<(v 2 ff) 2 > 



Measured 
value 



0.37 
160 
1.42 
0.43 

6.05 

0.02 
3.3 
0.49 ± 2% 



st ~ 0.01 

During an experimental run, the floating particles (50 
/im diameter and specific gravity of 0.25) are constantly 
seeded into the fluid from the tank floor, where they un- 
dergo turbulent mixing as they rise due to buoyancy and 
are uniformly dispersed by the time they rise to the sur- 
face. Once at the free-surface, their motion is constrained 
to the two-dimensional surface plane. Their motion is 
tracked with a high-speed camera (Phantom v. 5) situ- 
ated above the tank. The camera field-of-view is a square 
area of side length L = 9 cm. The constant particle injec- 
tion is necessary to replace floaters that stick to the tank 
walls. The sources and sinks at the surface fluctuate in 
both time and space, which can cause particles to leave 
the camera's field of view. 

Instantaneous velocity fields are measured using an in- 
house developed particle imaging velocimetry (PIV) pro- 
gram which processes the recorded images of the floaters. 
The constant injection of particles ensures that surface 
sources and sinks receive an adequate coverage of par- 
ticles on the surface. The local particle density at the 
surface determines the average spacing of the velocity 
vector fields produced by the PIV program. The result- 
ing velocity vectors are spaced (on average) by Sx = 2. 577 
over both sources and sinks , where 77 is the size of the 
smallest eddies in the inertial range [10| . 

Data were taken for several values of Re\ ~ 150 — 
170 with an average Re\ ~ 160. Turbulent parameters 
measured at the surface are listed in Table |TJ All of the 
statistics presented below were obtained by evolving ~ 
10 5 Lagrangian particles in each frame. 
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FIG. 3: A typical vorticity field displayed with the isolines. 
The image is of a square 9 cm by 9 cm. The thick solid line 
is the longest zero isoline in this particular field. 




200 400 600 800 



10° 10' 10 2 

FIG. 4: Main frame: a typical log-log plot of {(U(t + to) — 
[/(to)) 2 } showing its linearity with respect to t. Upper-left 
inset: a typical linear plot of {(U(t + to) — U(to)) 2 } showing 
the value of n as its slope. Lower-right inset: the distribution 

of k = {{8U{t))' 2 )/t 



IV. RESULTS AND DISCUSSION 

The vorticity field at each point, and at each instant of 
time r, is extracted from a measurement of the velocity 
field at equally spaced points separated by ~ 2.77/ = 0.54 
mm over the 9 cm x 9 cm field of view of the camera, 
giving a resolution of 160 x 160 points. The camera cap- 
tures the images at a rate of 133 Hz. The total number 
of frames that were analyzed is 833. 

In each frame there are roughly 1000 vorticity isolines 
and 400 zero- vorticity lines, see Fig. [3J The vorticity field 

is approximated using Stoke's Theorem lu — ^^ dl , where 



A is an area of size ~ (2.7rj) 2 . Although for a fractal 
object the concept of length is not clearly specified, the 
average length of the longest zero isoline of each frame is 
roughly 1000 steps, where the step size is 1.357y. 

From each frame, the longest zero and non-zero vor- 
ticity lines are extracted. An x-y axis is chosen for each 
isoline in each frame. The lowest point of the isoline is 
taken as the origin of the complex plane, x = y = be- 
cause gt(z) only concerns the upper complex plane, y > 0. 
Next, we calculate the driving function U(t) of each of 
the vorticity isolines using an algorithm implemented by 
T. Kennedy [ll[. 

Stated briefly, the algorithm goes by decomposing 
9s+t{z) = g~t ° 9s{z), where g a (z) maps the curve from 
to s back to the real axis and g~t(z) maps the im- 
age of the curve from s to t under g s (z) back to the 
real axis. By recursively decomposing g s (z), we get 

9t k = gk°9k-i°gk-2° — °52°<7ij0 = to < h < t2— < tk- 
Each cjk maps one small segment At of the curve to the 
real axis, where g~k has to satisfy Lowner equation (Eq. [2]) 
within each At. In principle, we take At — > but in the 
experiments the smallest At is the spatial interval of the 
two closest points with the same vorticity (At ~ 2.777). 
Each incremental map gk(z) is approximated by a Lau- 
rent series around z — > 00, gu{z) — z + ^ + Uk + O(-V)- 
This produces a sequence of discrete Uk that approxi- 
mates the true driving functions U(t). 

There are two ways of calculating ((SU(t)) 2 } (Eq. [1]). 
One way is to do self-averaging along t followed by an 
ensemble average over the 833 runs, another way is to 
calculate the ensemble average over t directly as done by 
BBCF. If U(t) is Brownian, as expected of SLE traces, 
the two methods should yield the same result for k and 
the exponent a. The resulting values of k and a from the 
two different methods will be denoted K t and at for the 
self-average and n e and a e for the ensemble average. 

The self-averaging is done by calculating ((U(t + t ) — 
[/(t )) 2 ), where (• • •) is an average over to for each 
fixed value of t. The exponent a in Eq. [T] is calcu- 
lated by taking the log of both sides of the equation, 
i.e. log((C/(t + t ) - ?7(to)) 2 ) = alog(t) +log(/t) for each 
curve and fitting it to a straight line (Fig. [4]). The slope 
of the line is the value of the exponent a for that partic- 
ular curve. The ensemble average over the 833 values of 
k and a are then taken to be the value for re t and at for 
the system. 

For SLE, as for Brownian motion, these two methods 
of averaging should produce the same results. Therefore, 
an easy way to distinguish vorticity data from effects of 
random noise follows from comparing results obtained 
by these two different averaging methods. To illustrate, 
we simulated SLE traces, and performed this particular 
test on those traces for a range of K-values from 1 to 
8 with increments of 0.5. In this simulation, there are 
1000 SLE traces for each value of n; each trace contains 
4000 points. We compare the results of ((SU(t)) 2 ) of the 
two averaging methods and find that they yield the same 
values of k within 2 %. 
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FIG. 5: Main frame: the probability distribution function of 
U(t)/(Kt) 1/2 for the zero-vorticity isolines for three different 
values of t. Lower-right inset: the log-linear plot of the PDF 
for 3 different values of t. 



TABLE II: Results of different tests of conformal invariance 
for zero isolines. For the non-zero isolines the results show 
that they are not conformally invariant 



Test 


t Avg 
(833 runs) 


Ensemble Avg 
(833 runs) 


SLE K 


Exponent a 


a t = 1.2 ± 0.1 


<z e =0.97 


a t =a e =l 


ft 


ftt =3.9 ± 1 


ft e =2.45 ± 0.2 


K± — Ke 


p(u(t)/Vrt) 




Not gaussian 

see Fig[5] 

and text 


<^)=0 


P{9) Eq. E] 
see dotted line 
in Fig. H 


No 


Yes 


see text 


D q independent 
of q, q=0 to 10 
see Fig. [7] 


see text 


No 


Yes 



We also extracted isolines from purely random fields. 
These fields were produced by generating uniformly dis- 
tributed vorticity amplitude at each site and have the 
same resolution as the measured velocity fields. The iso- 
lines of the random fields satisfy other tests of SLE (as 
described below) but fail the ensemble vs self average 
comparison test. Thus, comparing the results of the two 
averaging methods proves to be a very useful tool for 
differentiating noise effects from real SLE traces. 

Turning now to our experimental data, as seen in Table 
HIl the values produced by the two averaging procedures 
differ by approximately two standard deviations for both 
k and a. Since at and a e are measurably different from 
unity, the meaning of the k's in this case is ambiguous. 
We defer discussion of other parameters in the table. 

Two stringent tests of SLE, namely \a an d xl as de- 
scribed by Kennedy [l2j , were also performed on the iso- 
lines. The experimental data failed this test, and so did 



a simulated SLE at a similar resolution. Clearly a more 
refined resolution is required to produce a conclusive re- 
sult. 

The function U(t) is a random variable, and we have 
measured its probability distribution function (PDF, see 
Fig. [5]) to determine if it is gaussian for all values of t 
as required for a Brownian process. If so, the PDF of 
[/(^/(Kei) 1 / 2 should collapse onto a standard gaussian 
PDF with mean /i = and standard deviation a = 1 
(Fig. [5]). Here, only n e is meaningful, since the PDF of 
U(t) is produced by the values of U (t) for a fixed t in the 
ensemble. The inset of Fig. [5] shows that the PDF's of 
U(t) for three different values of t do not conform to a 
gaussian distribution. That is, all data points do not lie 
on an inverted V. For the non-zero isolines the PDF of 
Uit^Knet) 1 / 2 has non-zero mean and is strongly skewed, 
ruling out conformal invariance (mean and skewness of 
the zero isolines are 0.04 and -0.0006 respectively while 
for the non-zero isolines they are 0.6 and -0.2 respec- 
tively) . 

We apply another test to validate the Brownian char- 
acteristic of U(t) (Eq. Q}. If the vorticity isolines are 
characterized by Brownian U(t), they must be identi- 
fied by a particular probability distribution for keeping 
a point z = pe l9 to the right of each isoline. Here 9 is 
defined with respect to the x axis. It is given by .13;] 



P(0) 



i-Fi -; -;-cot 2 0j cot(0), 

(3) 

where T is the ordinary Gamma function with 
parameter, and 2-F1 is the gauss hypergeometric function. 

The probability distribution for the zero isolines is 
shown in Fig. [S] The solid lines represent the distri- 
butions based on the mean values of k; the dashed lines 
denote ± one standard deviation. The dots are the mea- 
sured distribution. The left panel of this figure shows 
better agreement with our measurements. Equivalently, 
the measured angular distribution fits better to the ex- 
pected distribution for n e . We are puzzled by this finding. 

Our last test to see whether the system exhibits confor- 
mal invariance is calculating the multifractal spectrum of 
the isolines (Fig. [7]), since conformal invariance requires 
scale invariance. The multifractal spectrum D q of the 
longest zero and non-zero vorticity isoline of each frame 
is then computed using 14] 



D„ = lim 

q^O q—1 



1 dlog(C g (r)) 



dlogr 



C q {r) 



1 N 

-y 



N-l 



N-l 



n 3 ) 



9-1 



(4) 



(5) 



Here N are the total number of points in the isoline, H 
is the heaviside step function, and r\j is the distance be- 
tween points i and j. This algorithm for determining the 
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FIG. 6: Probability that an isoline keeps a given z = pe to its right. The horizontal axis is 9 in radians with respect to the 
x axis. The probability distribution for zero isolines compared to the expected values for K e (Left) and Kt (Right). Solid lines 
are predicted distribution based on mean of k and dashed lines are based on the value of re± one standard deviation. 
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FIG. 7: Main frame: multifractal spectrum of the zero vortic- FIG. 8: Log- log plot of typical C q (r 



1/9-1 



versus r for q=2, 



ity isolines compared to the value of D q based on Kt- Lower- 
right inset: multifractal spectrum of the zero vorticity isolines 
compared to the value of D q based on K e 



5, and 10. The scaling region (straight portion) of C q {r) typ- 
ically extends 0.7 decades in r as depicted in the plot. From 
top to bottom, the curves are Cio, Cs, and C2 



spectrum of fractal dimensions is given by Hentschel and 
Procaccia 



15] 



To calculate D q , the log of the correlation sum (Eq. [5]) 
is plotted versus the log of r. The range of r over which 
the plot is a straight line is the scale-free (or scaling) 
region. The slope of the line d\og(C q (r))/d\og(r) is the 
value of D q . 

There are three lines in each graph in Fig. [7] They are 
the fractal dimensions given by D K = 1 + k < 8 [H, 
[l7j . The middle horizontal line in each graph represents 
D K using the measured mean value of k (taking a to 
be unity); the upper and lower dashed horizontal lines 
show D K using one standard deviation from the mean 
value of k. The dots represent the measured mean fractal 
dimension for each q. The vertical error bars show the 
uncertainty in the measured values of D q . The scaling 
region of C q (r), of which slope determines the value of 



D q , is shown in Fig. [3] The scaling region typically 
extends 7/10 of a decade in r for both q = 2 and q = 10. 

The multifractal spectrum of the zero isolines is then 
compared to that of the expected value of D K . Fig. [7] 
shows that the multifractal spectrum of the zero isolines 
conforms better to the value of D K given by n t . It is in- 
teresting to note that P(0) fits better to the distribution 
given by K e . Strictly speaking, a conformally invariant 
curve should have a constant multifractal spectrum which 
is independent of q. 

In light of the error bars in this figure, we are forced 
to conclude that departures from homogeneous fractal 
behavior are not clearly present. At the same time, the 
secular decrease in D q with increasing q, suggest that D q 
may not be a homogeneous fractal. 
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V. SUMMARY 

The goal of this experiment is to determine if the vor- 
ticity isolines of the compressible system of floaters meet 
all the imposed tests for SLE. The floaters exhibit a mea- 
surable departure from SLE for the zero vorticity isolines 
and a much larger departure for lines of nonzero vortic- 
ity. The departures are clearly evident in Fig. [H The 
interest of this study may be that the floaters display 



approximate SLE, as Fig. [4] and Table ILT1 show. 
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